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1. INTRODUCTION 

Electrical power generation map will change in the near future with lower environmental emissions. 
This will be done using renewable energy resources. Dependability on these resources is more and more 
increased whenever their operation problems more decreased. One of these problems is high variability 
which caused unstable systems. Solar energy is one of the major renewable energy resources used all over the 
world in electrical power generation. As penetration rate of solar power generation in electrical power 
systems is increasing with time, system stability is become important matter to be studied. Enhancement of 
renewable power operation stability will enhance the overall electrical power system stability with making 
the best use of these sustainable and clean resources. 

Electrical power generation from solar energy came in many forms; one of them is Photo Voltaic 
panels (PV). To maximize generated electrical power from a PV system, different maximum power point 
tracking algorithms (MPPTs) are used. These algorithms try to reach maximum power can be generated from 
the panels under valid solar irradiance and air temperature. Fractional open circuit voltage (FOC) [1], 
Incremental conductance(IC), Perturb and observe (P&O)[2] and Fractional short Circuit current (FSC) are 
examples of MPPT algorithms [3]. To enhance dynamic operation of different MPPT algorithms and giving 
better PV system dynamic response under heavy disturbances (variable irradiance and temperature), different 
types of controllers are used. A controller is a device used to give better stability and better system response. 
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Proportional-integral controller (PI) and Proportional-reference frame controller (PR) [4], Fractional 
order proportional-integral-derivative controller (FOPID) [5] Sliding mode control [6] and Fuzzy logic 
controller [7] are examples for controllers used to enhance operation of PV systems and MPPT. But these 
controllers need optimization to tune their gains. 

Optimization techniques are different algorithms used when minimization or maximization of an 
objective function is required. Different optimization techniques are used to tune controller gains for optimal 
maximum power point tracking (MPPT). As an example, in [8] Genetic algorithm and Particle swarm 
algorithm are used to track the MPP .Artificial Bee Colony algorithm is used in [9] for MPPT under partial 
shading conditions. In [10], authors showed a comprehensive review between different MPPT techniques for 
PV system output power maximization. In [11], authors made a practical performance evaluation for Fuzzy 
logic, P&O, Hill climbing and Incremental conductance algorithms for MPPT in a PV system. A complete 
comparison between IC and P&O algorithms is made under dynamic weather conditions in [12]. Authors of 
[13] made analysis of PV efficiency enhancement by using Incremental Conductance under non-linear 
loading conditions. A variable step modified fractional order IC technique is reported in [14] to maximize 
power tracking in fuel cells. Results showed better dynamic response than using default IC algorithm. 
Authors of [15] made simulation and experimental design of new advanced IC algorithm with variable step 
size for PV systems. Results showed better PV array output power response under sudden irradiance 
variations than using default IC algorithm. Enhancement of IC algorithm operation under partial shading is 
made using direct fuzzy controller in [16]. In [17], authors enhanced transient stability operation of 20 MW 
PV power plant by modified Proportional-Integral-IC algorithm controller and Fractional Order Proportional- 
Integral-IC algorithm controller for maximum power from the plant. 

In this paper, a 100 KW grid connected PV system in MATLAB SIMULINK is adopted to 
withstand sudden irradiance and air temperature variations due to partial shading on the PV panels. 
Incremental Conductance MPPT is the algorithm used to maximize system output power. New adaptive 
controller is used to enhance Incremental conductance operation to prevent MPPT failure. Two different 
Optimization techniques (Invasive weed optimization (IWO) and harmony search algorithm (HSA)) are used 
to tune the designing parameters of the new controller. Different optimization techniques results are 
compared. A Robustness test is made to check system ability to track maximum power point under six 
different random irradiance patterns. Finally, the PV system in [18] is constructed to compare the new 
optimized adaptive controller with the used P&O algorithm to track the maximum power point. 


2. ON-GRID PV SYSTEM MODEL IN MATLAB 
A 100 KW on grid PV detailed model in MATLAB SIMULINK is used. The model name is 
“power_PVarray_grid_det” and has connections as shown in Figure 1. The model main components are: 
a. PV array: To generate 100 KW electrical power 6 modules are connected in series to form a string and 
66 strings are connected in parallel. The module used name is “Sun power SPR-305-WHT” with 
electrical specifications as in Table 1. 


Table 1. DC Electrical Characteristics 


Quantity Value 
STC Power Rating 305W 
PTC Power Rating 280.6W 
STC Power per unit of area (187.0W/m’) 
Peak Efficiency 18.7% 
Power Tolerances -5%/+5% 
Number of Cells 96 
Nominal Voltage not applicable 
Imp 5.58A 
Vmp 54.7V 
Isc 5.96A 
Voc 64.2V 
NOCT 45°C 
Temp. Coefficient of Power -0.38%/K 
Temp. Coefficient of Voltage -0.177V/K 


b. Boost Converter: used to take PV array output DC voltage and current and change them to another DC 
voltage and current values to reach MPPT. Converter has a 500 volt max operating voltage and 5 KHZ 
switching frequency. Converter switching pulses came from converter controller circuit (Incremental 
conductance algorithm+ Integral regulator) 
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c. Inverter: takes converter output and starts to convert DC voltage and current into AC voltage and current 
using 3-level bridge circuit and a Pulse Width Modulation control circuit. 

d. Current Filter: To decrease total harmonic distortion in the AC current from the inverter and smooth 
current variations. 

e. Transformer: a “100 KVA-260V/25KV” transformer is used to step up voltage to the required grid 
voltage. 

The boost converter control circuit is modified from the default Incremental Conductance 
and Integral Regulator to new control circuit (Incremental conductancet+Adaptive Controller (Integral 
Regulatort+tSecond Order Amplifier)). The second order amplifier is connected as shown in 
Figure 2. The amplifier transfer function between input and output is in equation (1, 2).The amplifier 
designing variables for optimal adaptive controller tuning are charging time constant (Tp) and discharging 
time constant (T;). 
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Figure 1. On-grid-100 KW PV system model in MATLAB SIMULINK 
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Figure 2. Boost converter control circuit with incremental conductance MPPT and adaptive controller 


(output) _ 1 (1) 
(input) — Tp.T)s2+Tpst1 


Int J Elec & Comp Eng, Vol. 8, No. 5, October 2018 : 2650 - 2660 


Int J Elec & Comp Eng ISSN: 2088-8708 O 2653 


(output) _ 1 
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3. OPTIMIZATION TECHNIQUES 
Two Different optimization techniques are used in this paper for PV system better performances 
which are Invasive Weed Optimization (IWO) and Harmony Search Algorithm (HSA). 


3.1. Invasive Weed Optimization (IWO) 

Steps of Invasive Weed Optimization technique are [19]: 

a. Initializing a population: finite number of populations (seeds) is being selected with random positions. 
Reproduction: every seed grows to form a plant and produces newer number of seeds depending on its 
fitness value. 

b. Spatial dispersal: The produced new seeds are being randomly dispread over the search area and grow to 
new plants. Randomness and adaptation in the algorithm is provided in this part. The generated seeds are 
being randomly distributed over the d dimensional search space by normally distributed random numbers 
with mean equal to zero; but with variable variance. During simulation, a nonlinear alteration has shown 
satisfactory performance which is given in equation 3. Where iter, is the maximum number of 
iterations, iter is the present time iteration, oj¢¢, is the standard deviation at the present time step and n is 
the nonlinear modulation index. 

c. Competitive exclusion: In this process maximum number of plants is reached and then the process 
continues until maximum iteration is reached and hopefully the plant with best fitness closest to the 
optimal solution is obtained. 


_ (itermax—iter)” 
r (Ginitiat 


Oiter = z Osina) F Ofinal (3) 
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3.2. Harmony Search Algorithm (HSA) 
Harmony Search Algorithm to maximize certain objective function has the following steps [20]: 
1. Initialization 
a. Setting design variables array with dimension N: [X1, X2, ..., Xi ...-Xwy] 
b. Set population size and xX 1 is the population number "j? for the designing variable 
“i”; Where j=1, 2....., (population size). 
c. Initialize objective function f (X/) will be maximize 
d. Set each designing variable range X} < X; < X¥ and start with random point and dimension N. 


xÍ = X} + (XË — X}) + rand(1,N) (4) 


e. Set HMCR value which is a probability function takes value between 0 and 1. 

f. Select Pitch Adjustment Rate (PAR) with ranges from 0 to 1 (PAR=1 is used). 
g. Select the step size (b;) which is fixed or computed number (b; = 1 is used). 
h. Set iterations maximum number (used=25) 

2. Initiation of Harmony Memory (HM) 


Xi Xi m XG FR) 
2 2 2 2 
HM= Xi X2 wo XN (5) , F= f(x") (6) 
XHMS yHMS | yHMS f (XMS) 


3. Generation of New Harmony is based on the following criteria 


a. New Harmony is selected randomly; X’ = (Xj, X3,........Xy) with a (HMCR =1) probability 
b. Choosing probability (PAR = 1) doesn’t add any amount to New Harmony and New Harmony will be 


c. Worst value for the objective function is removed from harmony memory at end of iteration. 
4. Repeat till maximum iteration value is reached. 
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4. OBJECTIVE FUNCTION 

The Objective function will be maximized using both optimization techniques is the delivered 
energy (E) to grid form the 100 KW PV system. The delivered energy (E) is the integration of PV system 
delivered power to the grid (P) with time as in equation (7). Increasing the delivered energy means 
enhancment of incremental conductance algorithm to be stable under sudden irradiance and temperature 
changing patterns without failure to track the maximum power point (MPP). The designing variables for the 
second order amplifier to maximize delivered energy will be amplifier circuit charging time (T_h) and 
discharging time (T_l). Let O and M defined as in equation (8-9). Ranges used for O and M by the 
optimization techniques to maximize the objective function will be as in equation (10,11). The optimization 
algorithms will select optimal O and M values for objective function maximization. 


E = f; P. dt (KW.sec) (7) 
—_ ti -1 

0 = T (sec™>) (8) 
Zel -2 

M= FT (sec “) (9) 

0.1<0<50 (10) 

0.1 <M < 100 (11) 


5. METHODOLOGY 

Optimization techniques (IWO and HSA) general working methodology to select optimal charging 
time and discharging time values for optimal tuning of the adaptive controller to maximize the delivered 
energy (E) to the grid is shown in Figure 3. 
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Figure 3. Flow chart of methodology 


6. SIMULATION RESULTS AND COMPARISON 

The simulation results are divided into three stages. First one is PV system response under (variable 
irradiance pattern—constant PV cells temperature). The second one is PV system response under (variable 
irradiance pattern—variable PV cells temperature). A Robustness test is carried out to check the stability of 
the new optimized controller. Finally; the new optimized controller is used with PV system constructed in 
[18] and compares its performance with the used P&O one. 
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6.1. PV System Response Under (Variable Irradiance Pattern-Constant PV Cells Temperature) 

Irradiance pattern shown in Figure 4 is applied as an input to the PV array without taking cell 
temperature variations into consideration. Cell temperature is assumed constant as in Figure 5. The 
optimization process is made one time using invasive weed optimization algorithm and another time using 
harmony search algorithm to tune the gains of the optimized controller. As shown in Table 2, Invasive weed 
optimization algorithm gives slightly higher delivered energy (E) to the grid but taking longer optimization 
time. While harmony search algorithm has very lower optimization time. Figure 6 shows the PV system's 
delivered power to the grid after and before using optimized controller. Before using the optimized 
controller, power decreased to zero in the time intervals from (t=1.18 to t =1.25 sec and t=1.5 to t =1.52 sec) 
but by using both IWO-optimized controller and HSA-optimized controller the system continues tracking the 
MPP without falling to zero at these time intervals. 


Wr aan etim’) 


Figure 4. Solar irradiance pattern with time 
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Figure 6. Delivered power to the grid with and without optimized controller 
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Table 2. IWO and HSA Optimization Results for Constant Cell Temperature 
Optimization TWO HSA 
Charging time inverse value (0) (sec~*) 3.004 9.277 
Multiplied charging and discharging times 
inverse value(M) (sec™?) 


30.648 24.0713 


Objective Function( E)( KW.sec) 139.8616 139.5044 
No. of Iterations 25 25 
Time Taken(hours) 26.2 2.3 


6.2. PV System Response Under (Variable Irradiance Pattern-Variable PV Cells Temperature) 

The optimization process is repeated again but with taking cell temperature variations into 
consideration. The assumed temperature pattern is seen in Figure 7. The temperature pattern is constructed 
within a correlation coefficient range (between irradiance and air temperature variations) as in reference [21]. 
As shown in Table 3, IWO also gives slightly higher objective function value (E) but with longer 
optimization time.HSA also has very lower optimization time. The PV system's delivered power to the grid 
after and before using optimized controller is shown in Figure 8. 

The system before optimization failed to withstand sudden decrease in the irradiance and system 
output power dropped to zero in the time interval from (t=1.18 to t =1.71sec) but with I[WO-optimized 
controller the system continues MPPT. HSA-optimized controller also successes to track the MPP at that 
time interval. However IWO gives slightly delivered energy to the grid but HSA is more robust against 
temperature and irradiance variations. It's cleared from the previous results that without the optimized 
controller there is a MPPT failure. 


Table 3. Optimization Techniques Results for Variable Cell Temperature 
Optimization IWO HSA 
Charging time inverse value (0) (sec™*) 4.643 4.464 
Multiplied charging and discharging 


35.252 26.146 
times inverse value(M) (sec ~?) 
Objective Function( £)( KW.sec) 138.379 138.299 
No. of Iterations 25 25 
Time Taken(hours) 28 2.5 
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Figure 7. Variable cell temperature pattern Figure 8. Delivered power to the grid with and without 


optimized controller with variable cell temperature 


6.3. Robustness Test 

To check the robustness of the optimized controller six random irradiance and temperature patterns 
are used. Patterns are applied to the system without repeating the optimization process. From Figures 9 and 
10, [WO and HSA proved their reliability through keeping the system at the the MPP although the normal 
system without optimized controller failed in that. 
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Figure 9. Robustness test for three different irradiance patterns with corresponding temperature variation patterns 
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Figure 10. Another different three irradiance patterns with corresponding temperature patterns for more convenient 
robustness test 
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Figure 10. Another different three irradiance patterns with corresponding temperature patterns for more convenient 
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6.4. Comparison 

In this section, the more robust technique shown in the previous results (HSA-optimized controller) 
will be compared with the used P&O algorithm in reference [18]. The PV system is constructed as in 
reference [18] using MATLAB/SIMULINK with two PV arrays with about 200 KW rating. The gains were 
obtained in the first run optimization will be used without change to proof that the HSA-optimized controller 
is still reliable even with changing the PV array rating. The input irradiance for the first and the second PV 
array is shown in Figure 11. 

The temperature is assumed constant at 25C as in [18]. The power delivered to the grid obtained in 
reference [18] is shown in Figure 12 while P&O algorithm was used to track the MPP. Figure 13 shows the 
delivered power to the grid with the HSA-optimized controller to track the MPP. From Figures 12 and 13 the 
area under the curve for the delivered power of the HSA-optimized controller is higher than area under the 
curve of the P&O algorithm, so HSA-optimized controller gives higher delivered energy to the grid which 
means better average MPPT. Finally from the obtained results, robustness test and comparison, HSA- 


optimized controller proved its reliability under variable irradiance, variable temperature and variable PV 
rating. 
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Figure 11. Irradiance patterns for array one and two Figure 12. Delivered power to the grid using P&O 
algorithm 


Figure 13. Delivered power to the grid using the HSA-optimized controller 


7. CONCLUSION 

Heavy sudden disturbances such as sudden irradiance changes due to partial shading and 
corresponding temperature variations may result in MPPT failure and make PV systems unreliable. Under 
different random irradiance patterns the incremental conductance with enhancement of integral regulator 
failed to withstand these disturbances. Using the new optimized adaptive controller consists of integral 
regulator and second order amplifier results in no MPPT failure so; enhancement of incremental conductance 
operation and more PV system reliability is achieved. After making the optimization process to tune the 
adaptive controller gains better results are reached. The IWO-optimized controller gives slightly higher 
delivered energy to the grid than The HSA-optimized controller but with longer time for optimization. The 
HSA-optimized controller is more robust against irradiance and temperature variations. Finally from the 
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obtained results, robustness test and comparison HSA-optimized controller proved its reliability under 
variable irradiance, variable temperature and variable PV rating. 
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